Preprint typeset in JHEP style - HYPER VERSION 



o 
o 



m 



o 



o 



Z(3)-symmetric effective theory for SU{3) 
Yang-Mills theory at high temperature 



A. Vuorinen* and Laurence G. YafFe^^ 

Department of Physics, University of Washington, Seattle, WA 98195-1560 



(N 

Q I Abstract: A three-dimensional effective theory for high temperature SU{3) gauge 



theory, which maintains the Z{3) center symmetry of the full theory, is constructed. 



t:J- I Such a Z{3) invariant effective theory should be applicable to a wider temperature 



range than the usual effective theory, known as EQCD, which fails to respect the 



O I center symmetry. This center-symmetric effective theory can reproduce domain wall 

i-C I and phase transition properties that are not accessible in EQCD. After identifying a 

I I convenient class of Z{3) invariant effective theories, we constrain the coefficients of 

Q I the various terms in the Lagrangian using leading-order matching to EQCD at high 

'^ ' temperature, plus matching of domain wall properties in the full theory. We sketch the 

expected structure of the phase diagram of the effective theory and briefly discuss the 

rS ' prospects of numerical simulations and the addition of quarks. 
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1. Introduction 

The thermodynamics of high-temperature QCD with various numbers of massless 
quarks has been widely studied using the method of dimensional reduction, in which one 
exploits the decoupling of the non-static degrees of freedom to build an effective three- 
dimensional theory known as electrostatic QCD (EQCD) |l|. This approach has been 
used to derive the perturbative expansion of the QCD pressure through order g^lng, 
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and to compute various gluonic and mesonic correlation lengths 0, ||, ||. It has also led 
to improved understanding of the non-perturbative dynamics on the "ultra-soft" g^T 
energy scale . This approach is applicable at asymptotically high temperatures where 
the running coupling g'^{T) is small. But closer to the confinement /deconfinement tran- 
sition, at temperatures of a few times T^, even non-perturbative simulations of EQCD 
have failed to produce satisfactory results [^ |^. There have been attempts to build 
effective theories for the Wilson line near the phase transition f^, |], |lO|, 11 1, but none 



of these theories (which often come with very complicated Lagrangians) provide any 
connection to the perturbative regime at T 3> T^. 

Our goal is to formulate a simple three-dimensional effective theory for hot QCD 
that can reproduce equilibrium properties at temperatures ranging from asymptotically 
large down to Tc. For simphcity, most of our discussion will focus on the case of pure 
Yang- Mills theory, although the addition of dynamical quarks will be briefly discussed 
in the last Section. The construction of our theory relies on a separation between the 
inverse correlation length and the lowest non-zero Matsubara frequency, 27rT. In SU{3) 
Yang-Mills theory, there is a parametrically large separation between these scales at 
very high temperatures, and a modest separation near Tc.^ 

In pure SU{3) Yang-Mills theory, the confinement /deconfinement phase transition 
is a symmetry-breaking phase transition. The theory has a global Z{3) symmetry cor- 
responding to invariance of the action under gauge transformations which are periodic 
in the temporal direction only up to a twist belonging to the center of the gauge group. 
This Z{3) symmetry is spontaneously broken in the deconfined phase where the Wil- 



^It is natural to distinguish the correlation lengths of time-reflection even and odd operators. The 
Debye screening length m^^ is the longest correlation length of time- reflection odd operators ||l3, such 
as the imaginary part of the trace of the Wilson line. The longest correlation length of time-reflection 
even operators, which we will denote by /i~^, is longer than the Debye length. In SU{3) Yang- 
Mills theory at asymptotically high temperatures, m'o/{2TTT) approaches g{T)/{2Ti) while /i/(27rT) is 
approximately as 0.1 ^(T)^ [^. (This is an estimate based on the corresponding SU{2) result, plus 
the fact that /i should scale nearly linearly with Nc-) Correlators of purely magnetic operators fall off 
like e"'^'''' at large distance, while correlators of the real part of the trace of the Wilson line initially 
decrease like e~^™i^'''l but eventually fall more slowly as e"''''''. The correlator of the imaginary part of 
the trace of the Wilson line initially decreases like e~'^™°'^' but switches to e~'"°'''' at large distance. 
As the temperature decreases, the ratio ii/{2ttT), as deduced from the Wilson line correlator, rises 
to about 0.5 at T ~ 2Tc Q || and then falls to about 0.1 just above T^ ^. The ratio ■mr,/{2TTT) 
on the other hand rises to near unity at T ^ 2Tc [0. Because the confinement/deconfinement phase 
transition in SU{3) Yang-Mills theory is only weakly first order [p6|, this ratio should decrease, similar 
to the behavior of /i/(27rr), closer to the phase transition, but we are unaware of any numerical simu- 
lations confirming this. Near T^ correlations of purely magnetic operators appear to fall exponentially 
with a length scale near (6T)^^ but should eventually switch to the longer scale /Li~^ seen in the Wilson 
line correlator |Il4|. 
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son line operator (which is an order parameter for this symmetry) acquires a non-zero 
expectation value. Close to the phase transition region, fluctuations in the Wilson line 
become increasingly important and the Z{3) invariance of the Wilson line probability 
distribution is an essential feature of the theory. Even in full QCD with dynamical 
quarks it has been argued that the Z{3) symmetry, although explicitly broken by the 
presence of light quarks, may nevertheless play an important role in the dynamics near 

The standard high temperature effective theory, EQCD, cannot provide a smooth 
interpolation between asymptopia and T^.. This is evident even without considering 
numerical evidence, as the construction of EQCD violates one of the central tenets of 
effective field theory — it does not respect all the symmetries of the underlying theory. 
EQCD breaks the Z{3) center symmetry explicitly, even in the absence of quarks. Its 
leading order Lagrangian is obtained by expanding the one-loop effective potential of 
the Wilson line Ijl^ around one of its three degenerate minima. Effects of the explicit 



Z{3) breaking are clearly visible in the phase diagram of the effective theory 0.^ 

An effective theory that preserves the Z{3) symmetry can capture the dynamics 
of thermal fluctuations (above T^) which create bubbles sampling all three degenerate 
equilibrium states, and thus should have a phase structure similar to that of the full 
four-dimensional theory. In particular, its phase diagram should contain a quadruple 
point, where all three deconfined phases coexist with the confining phase, thereby 
extending the description of the theory into a temperature range inaccessible to EQCD. 
Our guide during the process of constructing a better effective theory for A^f = 
QCD is the requirement that the new theory provide a smooth link between asymp- 
totically high temperatures, where perturbation theory is valid, and Tc, where a Z{3) 
symmetry restoring phase transition occurs. At high temperatures, its predictions must 
reproduce those of EQCD (within its domain of applicability). But a Z{3) invariant 
effective theory should also reproduce the thermodynamics of domain walls which in- 
terpolate between different equilibrium states related by the Z{3) symmetry |T^. 



We begin in Section 2 with a brief review of the relevant aspects of the thermo- 
dynamics of SU{3) Yang-Mills theory, including dimensional reduction and the role of 
the Wilson line as an order parameter for the confinement /deconfinement phase tran- 
sition. Section 3 presents the general form of our effective theory. We perform lowest 
order matching to ordinary A^f = EQCD at high temperatures in Section 4, and also 
compute the semiclassical domain wall profile in our effective theory. The expected 



^Even though it has been argued that a "partial dynamical restoration" of the Z{3) symmetry takes 
place in EQCD |^, this theory does not have three distinct, degenerate, and physically equivalent 
equilibrium states which would signal a spontaneously broken Z(3) symmetry. These aspects will be 
discussed in more detail in Sec. 5. 
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phase diagram of our theory is discussed in Section 5, and compared to that of EQCD. 
Possible future directions, including the addition of quarks, are briefly discussed in the 
concluding Section 6. Some details of the analysis of Section 3 are relegated to the 
appendices. 

2. Thermodynamics of SU{3) Yang-Mills 

2.1 Z{3) center symmetry 

The functional integral representation for the partition function of SU{3) Yang-Mills 
theory involves an integral over all gauge fields which are periodic in (Euclidean) time 
with period f3. In addition to invariance under gauge transformations which are periodic 
in time, the action and the space of field configurations are invariant under gauge 
transformations g which are only periodic in time up to a "twist" belonging to the 
Z{3) center of the gauge group, 

A^{x)^s{x){A^{x)+id^)s{x)\ s{x)eSU{3), (2.1a) 

s{x + (3 it) = z s{x) , zeZ{3). (2.1b) 

Under such a transformation, the trace of the Wilson line (or Polyakov loop) 



tr n(x) = tr < P exp 



i I drAo(r, x) 





(2.2) 



transforms in the fundamental representation of Z{3), 

tr n(x) -^ ztr l](x). (2.3) 

The magnitude of the expectation value of trr2(x) may be interpreted physically as 
the exponential of the change in free energy due to the addition of an infinitely heavy 
fundamental representation test quark at position x, 

|(trfi(x))| =e-^^^«("). (2.4) 

This expectation value vanishes in the confining low temperature phase, where the 
free energy cost to introduce a single test quark is infinite, but the expectation value 
is non-zero in the deconfined high temperature phase, where the free energy cost of 
a test quark is finite. Hence, the Z{3) center symmetry is unbroken in the confining 
phase, but is spontaneously broken in the deconfined phase, and the expectation value 
of tr r2(x) is an order parameter for the Z{3) center symmetry.^ 

■^Unlike more typical symmetry breaking phase transitions, the low-temperature (confining) phase 
is the disordered phase, while the high-temperature (deconfined) phase is the ordered phase. 
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Symmetry considerations [jT9|, as well as lattice simulations |T^, show that four- 



dimensional SU(3) Yang-Mills theory belongs to the same universality class as the three- 
dimensional three-state Potts model, and undergoes a weak first order phase transition 
at the critical temperature. At the phase transition, the theory has a quadruple point 
with four co-existing equilibrium phases (the three spontaneously broken deconfined 
states, related by Z{3) transformations, plus the confining phase). No correlation 
lengths diverge as T ^ Tc. 

The transition from the deconfined plasma above Tc to the confining glueball phase 
below Tc may occur through two mechanisms: bubble nucleation or complete wetting. 
In the former, bubbles of the confining phase are formed via thermal fiuctuations inside 
the (supercooled) deconfined matter just below T = Tc, and as the temperature is 
lowered these bubbles grow, eventually filling all space. Alternatively, complete wetting 
takes place if one begins with domain walls separating deconfined phases with differing 
Z{3) "magnetization" {i.e., differing phases for {tiQ)). Such a domain wall may be 
viewed as a thin film of the confining phase residing between the two deconfined phases. 
As the temperature is lowered below Tc, this film widens and splits the "deconfined- 
deconfined" domain wall into two "confined-deconfined" domain walls which repel, 
leading to the expansion of the confining phase in the direction perpendicular to the 
original domain wall. 

Domain walls separating Z{3) rotated deconfined phases are topologically stable, 
non-perturbative objects which play an important role in the dynamics of quarkless 
QCD at temperatures just above Tc. At sufficiently high temperatures, their proper- 
ties may be computed using semi-classical techniques |]T8[. In hot QCD, the effective 



potential for the Wilson line only arises at one loop order |T^. As a result, the De- 
bye mass rriD (which is determined by the curvature of this potential at its minimum) 
is of order g{T)T, while the energy density at the top of the barrier separating Z{3) 
minima is 0{T'^) (not order T^/ g{TY, as one might have expected for a semiclassical 
object) . The width of a domain wall is determined by the Debye screening length. So 
high temperature domain walls have an 0{T^) energy density over a distance of order 
m^, yielding a tension which scales as T^/g(T). This parametric dependence must be 
reproduced by any effective theory that claims to capture the Z{3) symmetry breaking 
physics of hot Yang-Mills theory. 

2.2 Dimensional reduction and EQCD 

Due to the compact (Euclidean) temporal direction in a finite-temperature field theory, 
all four-dimensional fields may be decomposed into Fourier sums running over Matsub- 
ara frequencies. Each Fourier component acts as a three-dimensional field, with a mass 
(for bosonic fields) given by an integer multiple of 2ttT. At sufficiently high tempera- 
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tures, all non-static modes act like heavy fields and decouple from dynamics on length 
scales large compared to 1/T. The remaining static modes are the relevant degrees of 
freedom for physics on longer distance scales, yielding an effective three dimensional 
theory. To lowest order, the form of the effective theory may be found by integrating 
out the heavy fields explicitly [0. At higher orders, it is more efficient to start with the 



most general form of the Lagrangian and adjust its parameters by demanding that the 
effective theory reproduce a minimal set of physical quantities which are computable 
in both the effective and underlying theories |jl|, pO| . 

In the case of SU{3) Yang-Mills theory at high temperature, the degrees of free- 
dom of the effective theory correspond to the static components of the gluon field A^. 
Integrating out non-static fiuctuations generates an order gT mass for the temporal 
component Aq, due to Debye screening. In addition, self-interaction terms are gener- 
ated for the Aq field. The resulting Lagrangian of A^f = EQCD reads^ 

-Ceqcd = 93^ { \ tr 4 + tr [( A ^o)'] + m| tr (A^) + Ae tr (4) } + 5£e , (2.5) 

with g^, = g{T)\/T the three-dimensional gauge coupling and Di = di — i[Ai, •] an 
adjoint representation covariant derivative. The final term (5£e stands for higher di- 
mension operators whose effects, on length scales large compared to T~^, are sub leading 
and only affect the thermodynamics at order g'^ . The operators in (5£e will be irrelevant 
for our purposes. 

A straightforward way to obtain the form of the above Lagrangian, plus the leading- 
order values of its parameters, is to consider the one-loop effective potential for the 
Wilson line r2(x) = Ve^Jo drAoiT,:x.) j^ ^j^g ^^jj theory. The resulting potential jl^ has 



three degenerate minima related by transformations belonging to the Z{3) center of 
SU{3). Choosing to work around the minimum where the Wilson line expectation value 
is real and positive (and subtracting from the potential an Aq contribution originating 
from fluctuations in the static gauge fleld, which are not to be integrated out in the 
present case), one obtains Eq. ( ^.5| ) as the resulting effective potential for Aq together 
with the lowest-order values for the effective theory parameters, 

ml = g{TfT\ \^ = \g{Tf/T:\ (2.6) 

To lowest order, the mass parameter rriE is the same as the physical Debye mass ttid- 
This derivation also reveals the major shortcoming of the effective theory: when 
constructing the EQCD Lagrangian, one has chosen to work in the neighborhood of 



''This form of the effective theory is vahd for SU{2) and 5*^7(3). With four or more colors, there 
is an additional (trv4§) term. But for both SU{2) and SU{2>), the two fourth-order terms are not 
independent, as (trAg) — 2tr(74g). 
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Ao ^ 0, corresponding to a real, positive Wilson line, and thus loses the Z{3) in- 
variance of the original theory. The other Z{3) minima correspond to eigenvalues of 
Aq near ±|7rT, and are outside the domain of validity of the effective theory. At 
asymptotically high temperatures this is of little consequence for most physical quanti- 
ties, as the probability of thermal fluctuations crossing the barriers separating different 
minima is exponentially small [with an exponent scaling as l/g{T)]. But fluctuations 
between the different minima become crucial for the dynamics closer to T^, and hence a 
non-Z(3)-invariant effective theory cannot properly describe physics near the confine- 
ment/deconfinement phase transition. 



3. Z{3) invariant effective theory 

3.1 Degrees of freedom 

To build a three-dimensional Z{3) invariant effective theory, the minimal degrees of 
freedom are the spatial gauge field A(x) and the Wilson line f2(x). But because ^(x) 
is a unitary matrix, a theory with polynomial interactions for the Wilson line will not 
be perturbatively renormalizable. For an effective theory which is only intended to be 
valid below a UV cutoff, this is not a fundamental problem. However, it is a serious 
practical nuisance. It complicates perturbative matching to the underlying fundamental 
theory, as well as matching between lattice and continuum regulated versions of the 
effective theory. 

Therefore, we take as our goal the construction of an effective theory whose short 
distance behavior is as benign as possible. To do so, we will replace the unitary matrix 
r2(x) by an unconstrained 3x3 complex matrix 2(x). This should be viewed as 
analogous to the relation between non-linear and linear sigma models (or between an 
Ising model and a double- well 0*^ field theory). One may view the complex field 2(x) as 
the result of applying a block-spin renormalization group transformation which averages 
r2(x) over a small spatial region, producing a result which is no longer unitary. 

An arbitrary complex matrix may always be decomposed into a product of unitary 
and hermitian matrices, so we may write 2(x) = r2(x) H(x.) with -ff(x) = /J(x)^ One 
may regard the matrix H{x) as containing unphysical degrees of freedom (not present 
in the underlying Yang-Mills theory) which we have chosen to "integrate in" in order 
to construct an effective theory with good short distance behavior. We will arrange for 
these extra degrees of freedom to be "heavy" — to have masses large compared to the 
inverse correlation length. 

We require that the effective theory Lagrangian be invariant under three dimen- 
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sional SU{3) gauge transformations, 

Z(x) ^ s(x) Z(x) s(x)t , 
A(x)^s(x)(A(x) + 2V)s(x)t, 

for s(x) G SU{3), as well as global Z{3) phase rotations, 

Z(x) ^e2™/3Z(x), (3.2) 

for integer n. 

Instead of using the polar decomposition of Z{x.) mentioned above, it will be more 
convenient for our purposes to decompose this field into a traceless part and its trace 
(times a unit matrix). Therefore, we define 

L(x) = trZ(x), (3.3) 

and 

M(x) = Z(x)-|trZ(x)l, (3.4) 

so that 

Z(x) = M(x) + i L(x) 1 . (3.5) 

3.2 Effective Lagrangian 

The Lagrange density for the effective theory will have standard derivative terms plus 
a potential for the complex scalar Z, 

C = gf {itr^ + tr [D.Z^D^Z) + V{Z)] , (3.6) 

with Di = di — i[Ai, ■ ] and Fij = diAj — djAi — [Aj, Aj\.^ The Z{3) invariant potential 
V{Z) will consist of two pieces. 



with 



and 



V{Z) = V,{Z)+glVr{Z), (3.7) 

Vo{Z) = ci tr [Z^Z] + C2 (det[Z] + det [Z^] ) + Cg tr [{Z^ Zf] , (3.8) 

Vi{Z) = ciii[M^ M] +C2{ii[M^]+ii[{M^f]) + c^ii[{M^ Mf] . (3.9) 



^The freedom to make multiplicative rescalings of Z is used to fix the coefficient of tr (DiZ^ DiZ) 
in the Lagrangian (3^). One could also add a |V(trZ)p term, as this preserves all the required 
symmetries. But this additional term may be eliminated with a -E — » Z + a(trZ)l field redefinition. 



We will refer to Vq{Z) as the "hard potential" and Vi{Z) as the "soft potential". This 
terminology is appropriate at high temperatures where g'^/T is small, in which case 
the effective masses arising from Vq{Z) will be large compared to those arising from 
9lVi{Z). 

The hard potential Vq{Z) is composed of the three lowest- dimension operators 
which (in addition to being 2'(3) invariant) are invariant under the S'f/(3) x SU{'i) 
transformation 

Z(x)^AZ(x)5, (3.10) 

for arbitrary A,B^ SU{3). Note that the free tr (diZ'^diZ) kinetic term is also in- 
variant under this transformation. A short exercise, presented in Appendix 0, shows 
that Vq{Z) has extrema when Z = XQ for Q any element of SU{3) and A a real root 
of 2c3 A'^ + C2 A^ + ci A = 0. The coefficient C3 must be positive for stability. We will 
require that C2 < and Cg > 9ciC3. These conditions ensure that the global minimum 
corresponds to the largest positive real root which we will write, for later convenience, 
as A = I f with 



-C2 + a/c|^^^~8c^ 



Therefore, Vq{Z) is minimized when Z equals an arbitrary SU{3) matrix times v/3.^ 
The coupling constants in the soft potential Vi{Z) will be required to satisfy the 
inequalities £3 > and C1C3 > c^. As shown in Appendix 0, these conditions ensure that 
Vi{Z) is strictly positive for non-zero M. Hence, the soft potential Vi{Z) is minimized 
when M = 0, or equivalently when Z is proportional to the identity matrix. The overall 
factor of gl multiplying the soft potential is designed to mimic the fact, noted earlier, 
that the effective potential for the Wilson line in high temperature Yang-Mills theory 
only arises at one loop order. 

We have only included in the soft potential the three lowest dimension Z{3) in- 
variant terms which can be built just from the traceless field M. There are more Z{3) 
invariant terms, with the same dimensions, which we could have added — for example, 
terms including factors of the trace L. However, if Z is restricted to the manifold of 
minima of the hard potential Vq, and so is proportional to an SU{3) matrix, then there 
are only three independent operators which are gauge invariant, Z{3) invariant, and 
at most fourth order in the field. The three terms we have included in the potential 



^Note that sending the coefficients Ci to infinity, while holding their ratios fixed, would cause Vq{Z) 
to act like a delta-function constraint forcing the matrix Z to equal an SU{3) element (up to a scale 
factor). The resulting effective theory would be a gauged non-linear sigma model for the Wilson line 
r2(x). But this limit would return us to a non-renormalizable effective theory, which we wish to avoid. 
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(|3.9|) are a convenient choice for these three independent terms. As shown in the next 
section, these terms will be sufficient for matching with the EQCD LagrangianJ 

The soft potential Vi{Z) lifts the continuous degeneracy at the minimum of Vo{Z), 
leaving three discrete minima related by the Z{3) symmetry. If Z equals | times an 
element of the center of SU{3), then it will simultaneously minimize both the hard 
and soft potentials, and hence will trivially minimize the combined potential V{Z) = 
Vo{Z)+glVi{Z). 

The theory defined by Eqs. ( p.6|) -( p^ is a superrenormalizable three-dimensional 
field theory. The only renormalizations required (in addition to vacuum energy sub- 
traction) are one and two loop adjustments of the coefficients ci and ci of the quadratic 
terms. 

To carry out perturbative matching, as discussed in the next section, we will need 
the expansion of the Lagrangian ( p.6| ) about the minima of the potential. In the vicinity 
of the minimum at | e^'^*"/^, we may decompose Z as 

Z = e2™/3 ji^i +^3 [75 (0 + ^x)l + (h + ta)] } , (3.12) 

where and x ^i^re real numbers, while h and a are traceless Hermitian matrices. We 
will also rescale the three-dimensional gauge field A —> (73 A to switch to conventional 
perturbative normalization of the gauge field. Note that the fluctuation fields 0, Xy 
h and a, as well as the gauge field A, will now have canonical dimension 1/2, as 
appropriate for a three dimensional theory. Expressed in terms of these fluctuation 
fields, the effective theory Lagrangian has the form 

C = Knin + |tr4 + 1 [{d,<pf + ml<p'] + \ [{daf + mix'] 

+ ti[{Dihf + mlh''] +tr[(Aa)'] + I4„t(0, X, Z^, o) , (3-13) 

where Di = di — igslAi, ■ ] and Vmin is the value of Vo{Z)/g'^ at its minimum. 

The interaction terms in ^nt are suppressed by one or more powers of gs, and come 
from the soft potential plus the parts of the hard potential which involve more than 
two fluctuation flelds. These terms are displayed explicitly in Appendix ^. The leading 
order masses induced by the hard potential are 

m| = -2ci - |f C2 , m1 = -vc2, m\ = -2ci - ^v C2 ■, (3.14) 



with V given in Eq. ( 3.11|) . The field a stays massless at this lowest order. Note that 



m\ = m^ + ml, so this mass is not independent of the other two. The relations between 
^Had we chosen to work with a gauge group SU{4) or higher, we would need to include a quartic 



double trace term in Eq. (3.9) in order to reproduce the EQCD Lagrangian. This is unnecessary for 



SU{3), due to the identity mentioned in footnote^. 
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the couplings in Vq{Z) and the tree level masses (and v) may be inverted to give 

Ci = |(m^ - 3ml) , C2 = -m^/t; , Cg = |(mj + 3m^)/f ^ . (3.15) 

The tree-level heavy masses m^ and m^, plus the value of v, should be regarded as 
arbitrary parameters for now, subject to the condition that m<^ > ^fri^. This is equiv- 
alent to the inequality Cg > 9ciC3 which ensures that the global minima of the hard 
potential Vq{Z) lie at |f e^'^*"'/^. The value of the potential (divided by g^) at these 
minima is Vmin = -f^ (9"^0 - ^\)/9l- 

4. Parameter matching at high temperature 

4.1 Reduction to EQCD 

We now want to deduce the constraints on the coefficients Cj in the soft potential 
Vi {Z) which must be satisfied if the effective theory is to reproduce the long distance 
dynamics of SU{'i) Yang- Mills theory at asymptotically high temperature. In this 
regime, fluctuations on the scale of the temperature T [as well as on the Debye scale 
g{T)T] are weakly coupled, so perturbative matching techniques may be used. The 
most efficient way to carry this out is to start with our 2'(3) invariant effective theory 
and integrate out the heavy fluctuation flelds 0, x ^i-nd h. This will produce a non-Z(3) 
invariant effective theory only involving the light field a (and the three dimensional 
gauge field), which may be directly compared to EQCD. 

To lowest non-trivial order, the result of integrating out the heavy fields is a "light" 
effective theory with Lagrange density^ 

Aight = I tr i^5 + tr[(A af + m^ a^ + A a^] . (4.1) 

Except for the trivial rescaling A ^ (73 A of the gauge fleld, and renaming Aq -^ g^ a, 
this has the same form, and identical derivative terms, as the EQCD Lagrangian ( p.5[) . 
The effective light mass m^, and the quartic coefficient A, can have contributions 
arising solely from the hard potential Vq{Z), as well as having contributions which 
involve the soft potential Vi{Z). The lowest order contributions generated by the hard 
potential are illustrated in Fig. |I|. Using the coefficients of the various terms in the 



^We omitted the unit operator from the original effective theory Lagrangian (3^), as weU as from 



(4.1), for simphcity of presentation. The coefficient of the unit operator gives the short distance 
contribution to the free energy density (divided by T). To match this contribution to the free energy, 
the coefficient of the unit operator in our effective theory must differ from the corresponding coefficient 
in EQCD by — Vmin + Y57 ("i^ + m^ + 8m^), up to higher order corrections suppressed by powers of 
(7I divided by a heavy mass. 
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a) 1« ^ X . 
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Figure 1: (a) Tree-level diagrams originating from the hard potential Vq which contribute 
to the quartic interaction vertex for the light field a. (b) One-loop diagrams produced by Vq 
which contribute to the mass of the field a. Solid, dotted, wavy and curly lines denote the 
propagators of the (p, x-, ^i and a fields, respectively. 

expansion of Vq about its minima, as given in Eq. (p.3| ), one may easily evaluate these 
diagrams explicitly. However, this is not necessary — the diagrams shown in Fig. |1|, 
at zero external momentum, sum to zero. This is a consequence of the SU{'i) x SU{?)) 
invariance ( |3.10| ) which the effective theory would have in the absence of the soft poten- 



tial Vi{Z) and the static gauge field. An infinitesimal SU{'i) x SU{'i) transformation 
acts as a shift on the light field a. Integrating out the heavy fields, without the soft po- 
tential, must produce a result which respects this invariance, and thus can only involve 
derivative interactions for the light field a. 

Consequently, the light mass m^ and the quartic coupling A in the Lagrangian 
(|4.1| ) can arise only from contributions involving the non-S'f/(3) x SU{?)) invariant soft 
potential Vi{Z). To lowest order, one may simply read off their values from the soft 
terms in the shifted potential ( B.3 ), giving 



^a = cig-i, X = C3g^. (4.2) 

Demanding that these coefficients coincide with the corresponding terms in the EQCD 
Lagrangian (OF) (after identifying Aq — > g^ a) fixes 



c, = T + 0{gl) , (4.3) 

and 

while C2 is left undetermined. 

Integrating out the heavy fields beyond leading order will, of course, produce cor- 
rections to the coefficients (^4.2|) , suppressed by additional powers of gl divided by a 
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heavy mass, as well as induce additional higher dimension operators in the light effec- 
tive Lagrangian ( |4.1| ). The first new operators, involving six powers of the light field 
a, will have a coefficient of order (s'l/mheavy)''- 

4.2 Domain wall properties 

Having ensured that our theory reproduces the perturbative physics of SU{3) Yang- 
Mills theory at very high temperature, we are left with two undetermined parameters, 
V and C2 (in addition to the deliberately introduced heavy masses m^ and rriy^). As we 
will show, these two parameters may be fixed by demanding that our effective theory 
reproduce the tension and width of domain walls separating different Z{3) minima. 

Consider a domain wall lying in the x-y plane, which interpolates between the state 
with {ti Z) = V at z = — oo and {ti Z) = e^'^'/^f at z = +00. To find the domain wall 
profile and evaluate its tension, we need to minimize the free energy for configurations 
of the scalar field Z{z) which depend on a single coordinate z and satisfy the boundary 
conditions 



L{z = —00) = V , 


L{z = 00) = e^"^' 


''v, 


(4.5a) 


U{z = -cx)) = 0, 


M{z = 00) = 0. 




(4.5b) 



To evaluate the (leading order) domain wall tension in the high temperature or weak 
coupling limit it will be necessary to include: 

1. The contribution of the hard potential, ^ J dz Vo{Z{z)). 

2. The gradient energy, ^ J dztrlDzZ^z)]"^. 

3. The contribution of the soft potential, T J dz Vi{Z{z)). 

4. The one-loop contribution from fluctuations in the gauge and scalar fields. 

Our task is greatly simplified by the hierarchical structure of the potential. Because 
of their overall factors of g^^, one might expect the first two contributions to be most 
important. However, a translationally non-invariant domain wall configuration can 
minimize Vo{Z), everywhere in space, provided Z{z) equals v/3 times some (spatially 
varying) SU{3) matrix. And the gradient energy will be comparable to the energy due 
to the soft potential, and not parametrically larger, if the width of the domain wall 
is inversely proportional to g^. Because the soft potential is multiplied by an explicit 
factor of g^ in the Lagrangian ( |3.6|) , its contribution is comparable to the effects due 
to fluctuations in the fields. Hence, the one-loop contribution cannot be neglected.^ 

^This parallels the computation of the leading order domain wall tension in the full theory ||18[. 
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Using gauge invariance, we may diagonalize Z{z) and write it as 

Z{z)=dmg{Mz)J,{z),fs{z)), (4.6) 

with 

/^(^) = |e¥H^)+3/3Wl, (4.7a) 

/2(^) ^v^'-r[aiz)-3m]^ (4.7b) 



■-47ri 



/3(^) = |e— "(^^ (4.7c) 

Here a and f3 are real functions of z that satisfy the boundary conditions 

a{z = — oo) = 0, a{z = oo) = 1, (4.8a) 

/5(^ = -oo) = 0, /3(z = cx)) = 0. (4.8b) 

Note that shifting {a,P) by any integer multiple of (|,±|) leaves 2{z) invariant (up 
to a permutation of its diagonal elements). 

It will be convenient to introduce a dimensionless rescaled coordinate z and Wilson 
line expectation value v via 

z = gsVfz, v = ^. (4.9) 

The gradient and soft potential contributions to the domain wall free energy (per unit 
area) then become 

/oo 
dz{ia'y + 3if3'y + U^ia,P)}, (4.10) 
-oo 

with a' = da/dz, etc., and 

Ui{a,(3) = — ^ js - cos(47r/5) - 2cos(27ra) cos(27r/3)| 

+ ^^ (s + cos(47ra) - 3 cos(47r/5) + 2 cos(27ra) rcos(67r/3) - 3 cos(27r/3)l | 

+ (^^) |9 + cos(87r/?) + 2cos(47r/3)[cos(47ra)-2] - 8 cos(27ra) cos(27r/?)|. 

(4.11) 

Fluctuations of the fields in the effective theory, at one-loop order, produce func- 
tional determinants depending on the background domain wall profile 2[z). The width 
of the domain wall will be parametrically large compared to 1/v, and therefore one may 
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use a gradient expansion in the evaluation of the resulting functional determinants. It 
is sufficient to keep just the first term, leading to^° 






\ndet{-d^ + [Z\ [Z, ■]]) /(Area) 

^ (1/1-/21^+1/2-/3^ + 1/3-/1^ 

/oo 
dzU2ia,/3), (4.12) 

■00 



with 



t/2(a,/3) = --^||sin(27r/?)|' + |sin(7r(«-/?))|'+|sin(7r(a+/?))|H. (4.13) 

Combining the three contributions gives 

Fdw[a, P] = -Fgrad + -^soft + -^fiuc 

/oo 
dz [(«')' + 3(/3')' + U^ia, /3) + f/2(a, /?)] . (4.14) 
-00 

This is the excess free energy density due to a domain wall (relative to the spatially 
constant equilibrium free energy density). Note that all three contributions are para- 
metrically the same size, as anticipated. Minimizing this functional, subject to the 



boundary conditions (O), will yield the leading order domain wall tension. 

The combined potential Ui{a,P) + U2{oi,[3) is plotted in Fig. 0. From the general 
form of the potential it is easy to see that for positive C2, and all values of a, there 
is a local minimum in the /^-direction at /5 = 0. This suggests, and further inspection 
confirms, that one may simply set /3 = in the search for the minimum. This reduces 
the computation to the minimization of 

/oo 
d-z[{a'f + U{a)l (4.15) 

-00 

^°The evaluation of this one- loop contribution is similar to the calculation in Ref. |Q. One may 
neglect the soft potential Vi{Z). Transverse fluctuations of A acquire an effective mass due to the 
background value of the scalar field, and yield the above functional determinant. If one uses the gauge 
fixing term tr (diAi)"^ and lets Z{x) — (1 + 93 r](x)) Z, with Z the background value of the scalar field, 
then the longitudinal part of A mixes with the traceless anti-Hermitian part of 77 (which receives no 
tree-level mass from the hard potential Vq). The resulting functional determinant is independent of 
Z, and equals one with dimensional regularization. Fluctuations of the trace and traceless Hermitian 
parts of 77 are also independent of Z, and produce an irrelevant constant (given in footnote @). 
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Figure 2: The potential Ui{a,(3) + U2{ct,l3) plotted for the optimal values of v and C2 
given in Eq. ( [4.231) . In the contour plot on the left, minima correspond to the dark regions. 
A fundamental domain, representing all configurations of Z which are inequivalent under 
conjugation by SU{3), is given by the triangle with vertices at (a, /3) = (0,0), (1,0) and 
(1/2, 1/2) — which represent the three distinct minima related by Z{3) transformations. The 
3d plot on the right shows the region between the two minima at {a, (3) = (0,0) and (1,0). 



with 



U{a) = — sin^(7rQ;) — 



2v 



sin (vra) + 



- sm vra 



(4.16) 



Solutions to the equation of motion produced by varying the functional ([4.15|) 
conserve the "energy" E = (S^a) — U{a), and the domain wall boundary conditions 
(|4.8|) imply that E = 0. The free energy functional (|4.15| ) is invariant under the 
transformation a{z) — » l — a{—z), which is a combination of parity, charge conjugation, 
and a Z[3) symmetry transformation. The minimal energy domain wall solution should 
have the same invariance, implying that a{z=0) = |. Hence, the appropriate solution 
to the zero energy equation of motion d^a = [/(o;)^/^ is given by 



a{z) 



daU{a)~^'^ = z, 



1/2 



while the resulting domain wall tension is 

a = Fdw[a,0] = g^^ (AnvTf {^Vff I da U{a] 



1/2 



(4.17) 



(4.18) 
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Figure 3: The function a{z) plotted for the optimal values of v and C2 given in Eq. ( [4.23 ). 



We will define the width Az of the domain wall as the ratio of the first moment of the 
domain wall free energy density divided by the tension/^ or 



Az 



(47rt;)2(iT)3 



^1^ 



dz \z\ U{a{z)) 



2{A7rvf arf '•I 



^1^ 



dai [/{aiy'"^ / da2 f/(«2) 

1/2 Jl/2 



-1/2 



(4.19) 



None of the integrals ( 4.17 )-( 4.19 ) are analytically tractable but evaluating them 
numerically, or integrating the equation of motion dza = U{aY^'^, for particular values 
of C2 and v is straightforward. The resulting tension and width may then be compared 
to the corresponding results in the underlying Yang-Mills theory. The (leading order) 
domain wall free energy density in high temperature SU{3) Yang-Mills theory is [J18l 

(4.20) 



•^ym(2;) 
The resulting domain wall tension is 



3cosh*(z(z)/2) 



ctym = / dzTYM(z) 



iTC^ 



J3 



9 9{T) 



and the (first moment) domain wall width equals 



A2;ym = o-vM / dz\z\J='^ 



"YMl^J 



ln(4) 



9iT)T 



(4.21) 



(4.22) 



^^Alternative definitions of the domain wall width could, of course, be chosen, but make little 
difference. For example, defining the width from the second moment of the domain wall free energy 



density changes the resulting optimal parameters (4.23) by about 0.1% for v and 2.6% for C2 
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Figure 4: Left: The domain wall free energy densities J^iz), divided by ivr^T^, in the effective 
theory and underlying Yang-Mills theory, plotted versus z = g(T)Tz. The two curves are 
nearly indistinguishable, differing by less than 1.5% at the center of the domain wall. Right: 
The relative difference of the two curves, (.Feff(z) — .Fym(^))/-^ym(^)- 

A straightforward calculation shows that the effective theory reproduces the domain 
wall tension and width of the underlying Yang-Mills theory when 



3.005868 , 



C2 = 0.118914, 



(4.23) 



[Note that this value for v, almost equal to 3, makes the minima of the potential lie 
almost exactly at T times an element of 2^(3).] In Fig. |^ we plot a{z) for these values of 
the parameters, and in Fig. ^ we compare the resulting domain wall free energy density 
in the effective theory, given by 



(4.24) 



J',^{z) = i4nvTy{lYU{a{z{z))), 
with the full theory result ( [4.201 ). The agreement is rather remarkable. 

4.3 Discussion 

Perturbative matching with EQCD, plus comparison of domain wall properties, has de- 
termined (to leading order) all parameters of the effective theory [leading to Eqs. ( [4. 3D , 
(|4.4| ), and ( t4.23| )] except the two heavy masses m^ and m^. These heavy masses are 
arbitrary, as long as they are parametrically of order T (or v), and satisfy m,^ > ^"rn^. 
It may seem surprising that domain wall properties can be used to fix parameters 
of the effective theory which were left undetermined by the perturbative matching to 
EQCD. Correct matching with EQCD should imply that our effective theory reproduces 
all physics of hot Yang-Mills theory on energy scales small compared to T. Fields in 
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the domain wall only vary on the Debye screening length scale [g{T)T]^^ which, for 
asymptotically high temperatures, is parametrically large compared to T^^. So why 
didn't matching to EQCD yield an effective theory that automatically reproduces the 
correct domain wall structure, without any further tuning of parameters? 

The resolution of this puzzle involves the realization that the inverse domain wall 
width is not the only relevant energy scale in a domain wall. As seen in the fluctuation 
contribution ( [4 .121 ) to the tension, near the center of the domain wall off-diagonal 
components of the static gauge field receive effective masses which are order v in the 
effective theory, and order T in the underlying Yang-Mills theory. ^^ Hence, the 0{T) 
value of V, which determines the magnitude of the scalar field in our effective theory, 
is a physically relevant scale in the domain wall. Equivalently, the excess free energy 
density ([4.24|) near the center of the domain wall is 0{T^), and thus comparable to the 
OiT"^) equilibrium free energy density due to fluctuations on the scale of T. 

The symmetry structure of the effective theory guarantees that it will have domain 
walls interpolating between different Z{3) minima, and the matching to EQCD guaran- 
tees that these solutions will have the correct exponential decay, proportional to e'^^kl^ 
far from the domain wall. However, the precise value of the domain wall tension, and 
its width, are sensitive to physics on the scale of T. This is why reproducing these ob- 
servables in the effective theory requires tuning of parameters beyond the perturbative 
matching to EQCD. If one were to demand that the effective theory exactly reproduce 
all moments of the domain wall energy density (|4.20|) , then it would be necessary to 



add an infinite number of additional higher dimension operators to the effective theory, 
with appropriately tuned coefficients. But from Fig. ^, it is evident that near perfect 
agreement is achieved, without introducing any higher dimension operators, just by 
suitably adjusting those low energy parameters {v and £2) of our superrenormalizable 
effective theory which were left undetermined by the matching to EQCD. 

To conclude this section, let us briefly review the different sources of corrections to 
the Lagrangian and parameters of the effective theory. These will become important if 
one aims at next-to-leading-order (NLO) accuracy in the determination of various phys- 
ical quantities through the effective theory, such as bulk equilibrium thermodynamics 
or the Z{3) domain wall tension and width. 



^^An alternative way to see this is to consider the electric field contribution to the Yang-Mills 
action, tr {DoAi)"^ — tr {d^Ai — ig{Ao, Ai])'^. A background value of Aq effectively shifts the Matsubara 
frequencies from 2TinT to 2nnT minus the eigenvalues oi gAo (in the adjoint representation). Near the 
center of a domain wall, the eigenvalues of Aq (in a gauge where Aq is static) are of order TrT/g — since 
this corresponds to an order one change in the phase of the Wilson line fi ^ qWI^Ao ^ Consequently, 
near the center of a domain wall there is no separation of scales between the effective frequencies of 
the static and non-static modes. 
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Integrating out the heavy fields, with vertices originating from the hard potential 
Vo{Z), will produce corrections to the scalar and gauge field kinetic terms. At first non- 
trivial order, the former involves determining the leading momentum dependence of the 
two-point graphs of Fig. |I|b, and will yield an order g"^ correction to the normalization 
of the scalar kinetic term in the Lagrangian of Eq. ( ^.1| ).^^ These corrections amount 
to multiplicative wave function renormalizations for the fields in question. 

Integrating out heavy fields at NLO will also generate corrections to the values of 
the parameters in the light theory (|4.1| ), whose matching to EQCD will require 0{g'^) 
relative adjustments to the couplings in the soft potential Vi{Z). For the parameters 
rria and A, one-loop contributions involving one soft vertex will produce 0{g'^) and 
0{g^) corrections, respectively. For bulk thermodynamics, the NLO correction to m^ 
is more important than the leading order value of A. 

Because the SU{3) x SU{3) invariance of Vq is not respected by the soft potential 
Vi, or by the coupling of the scalar fields to the gauge field, integrating out the heavy 
fields may also produce entirely new operators involving A and a of dimension three 
or higher. The first such new operators to appear include tra^ and tr (a^(Dja)^), and 
are suppressed by g^}"^ 

Ensuring that the domain wall properties of our theory agree with corresponding 
results in full Yang-Mills theory beyond leading order would, in addition to the above, 
require extending the derivative expansion of the one-loop fiuctuation contribution Fguc 
to next-to-leading order, as well as evaluating the two-loop effective potential for Z. 
The analogous NLO calculation in the full Yang- Mills theory is discussed in more detail 



in Ref. [21 1. The required analysis in the effective theory should be very similar. 



5. Phase diagrams 

5.1 EQCD 

The lack of Z{3) center symmetry in EQCD has been a primary motivation for the 
work presented in this paper. The consequences of this explicit breaking of the Z{3) 
symmetry are particularly apparent if one studies the phase diagram of EQCD for arbi- 
trary values of its parameters, as done in Refs. 0, ^ . One can form two dimensionless 

^•^Terms involving powers of the external momentum, which correspond to higher derivative terms 
in the Lagrangian, will be proportional to inverse powers of the heavy masses and are thus suppressed 
by additional factors of g when evaluated on momentum scales of order gT which are relevant in the 
effective theory. 

^^For tra^, one factor of g| comes from the necessary inclusion of one soft vertex from Vi. For 
tr (a^(_Dia)^), the multiplicative coupling constant is of order g^, but the covariant derivatives con- 
tribute additional factors of 173 — either explicitly or via spatial derivatives applied to gT scale physics. 
The same is true for other higher dimension operators involving derivatives. 
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Figure 5: The phase diagram of EQCD (from Ref. [Q) parametrized by x = 2'^E and 
y = m|/g'|. The phase with spontaneous symmetry breaking of the Aq -^ —Aq symmetry 
hes in the lower left region of the diagram; the upper region is the unbroken symmetry phase. 



The dotted "4d matching line" is the slice of the phase diagram which satisfies relation (5.1) 
giving the result of next-to-leading order matching with hot SU{3) Yang-Mills theory. 



ratios from the parameters (f^f, m|, and Ae) of EQCD, conventionally chosen to be^^ 
X = |Ae and y = m\/ g^. Parameter values which result from matching to hot S'f/(3) 
Yang-Mills theory correspond to the curve 



xy 



9 



87r2 



IGtt^ 



x + 0{x^ 



(5.1) 



where the next-to-leading order correction, proportional to x, has been included. In- 
creasing temperature corresponds to decreasing values of x (proportional to g(TY). 

For sufficiently small values of m'^ (or y) one finds that the Aq —>■ —Ao symmetry, 
which is a remnant of time refiection symmetry in the original 4-d theory, is sponta- 
neously broken. Within the broken phase, there are two degenerate equilibrium states, 
with opposite (and non-zero) expectation values for trAj]. The phase transition line 
separating the symmetric phase (with vanishing (tr Aq)) from this symmetry broken 



^^Our definition of Ae differs from that of Ref. [^ by a factor of 2/g. 
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phase is first order for small values of x and second order for large values, with a tri- 
er itical point at X ~ 0.26. This is illustrated in Fig. ^. On the first order part of the 
transition line there are three co-existing extremal equilibrium states, distinguished by 
positive, negative, and zero values of the order parameter (tr Aq). 

Near the first order transition line, these three equilibrium states with differing 
values of (tr A^) may be regarded as deformed remnants of the three equilibrium states 
of hot SU{'i) gauge theory (which are related by the center symmetry). But in EQCD, 
there is no three-fold symmetry relating these states, only a Z{2) symmetry inter- 
changing the two symmetry-broken states. The states with non-zero (trAj]) have, 
for example, energy densities and correlation lengths which differ from the symmetric 
{IyAD = state |[. 

Only the symmetric phase corresponds to the physics of the underlying hot Yang- 
Mills theory. However, as indicated in Fig. |^, the dimensional reduction line (|5.1j ) lies 
inside the broken phase! This line intersects the phase transition curve at x = (or 
T = oo), but for any non-zero value of x (or finite temperature T), matching with 
hot Yang-Mills theory requires studying the symmetric phase on the wrong side of the 
phase transition — that is, studying a metastable "supercooled" phase. For sufficiently 
small X, this is feasible. The lifetime of the metastable phase is so long that transitions 
to the stable symmetry broken phase are never seen in practical numerical simulations. 
But this is surely a less-than-satisfying feature of EQCD. 

As one increases x (or lowers the temperature) in EQCD, there is no region where 
a new "confining" phase appears, and no quadruple point where four equilibrium states 
coexist, analogous to the transition temperature Tc in real SU{?)) Yang-Mills theory, 
where the three spontaneously broken 2'(3) deconfined states and the ^(3) symmetric 
confining state all coexist. 

5.2 ^(3) invariant effective theory 

Viewed as a superrenormalizable three dimensional field theory with no constraints on 
its coupling constants (besides stability), our 2'(3) invariant effective theory (|3.6| - p.9|) 
depends on a total of seven parameters, which may be taken to be one overall scale f , 
the ratios of the heavy masses to this scale (m?/f^ and m^/f^), a dimensionless gauge 
coupling g1/v, and three dimensionless "soft" couplings (ci/f, 62, and Caf). Exploring 
the resulting phase diagram via numerical simulations, as all six dimensionless ratios 
are varied, is obviously impractical. However, it is inevitable that the parameter space 
of this theory contains (i) a region in which the Z{?)) symmetry is spontaneously broken, 
and (ii) a region in which the 2'(3) symmetry is unbroken. ^^ 
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Gauge symmetries are, of course, never truly spontaneously broken. In addition to the global 
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There are two ways in which one can easily imagine changing the parameters of 
the theory so as to move from the symmetry broken (deconfined, high temperature) 
phase to the unbroken (confining, low temperature) phase. If the heavy masses are 
varied and m^ becomes less than ^rn^, then the global minimum of the potential Vo{Z) 
switches from Z = vl (times a cube root of unity) to the symmetric point Z = 0.^^ 
If the theory is kept weakly coupled, so that g'^/v <^ 1, then this tree- level (or mean 
field) analysis is reliable, and the resulting transition will generically be a strong first 
order transition. 

Alternatively, if m^ remains greater than |m^ but the theory becomes sufficiently 
strongly coupled, so that g^/v is no longer small, then this should drive a fluctuation in- 
duced phase transition in which the entropy-driven effects of fluctuations overwhelm the 
symmetry-breaking influence of the potential. This is analogous to symmetry restora- 
tion in non- linear sigma models. 

The confining/ deconfining phase transition in SU{3) Yang-Mills theory is known 
to be a relatively weak first order transition, and near the transition fiuctuations of the 
Wilson line, or the static gauge field, are not weakly coupled. Consequently, we expect 
that the second regime of fiuctuation induced phase transitions is the more appropriate 
model for SU{3) gauge theory near its confinement transition. 

An interesting two-dimensional slice of the parameter space of our effective theory 
which would probe both of the above regimes, as well as the portion of parameter space 
which corresponds to asymptotically high temperature Yang-Mills theory, is given by 
varying g^/v and rn'i/v'^ while fixing the soft couplings at the values determined from 
the high temperature matching in the previous section, namely 

ci/t; = 0.332683 , £2 = 0.118914, 53^7 = 0.228419, (5.2) 

and choosing m'^/v^ to equal, say, unity. Within this two-dimensional slice, there will be 
a region of unbroken Z{3) symmetry, a region of spontaneously broken Z{3) symmetry, 
and a phase transition line separating these regions. The entire phase transition line is 
expected to be first order (for the same reasons discussed in Ref. [0), and everywhere 

Z{3) symmetry, our effective theory is also invariant under Z — > Z* (for real values of the couplings). 
Hence, the global symmetry group is really Z{3) x Z{2) = Z{6). The Z{3) broken phase is a phase 
in which this Z{6) global symmetry spontaneously breaks to a Z{2) subgroup. For real values of the 
couplings, we do not believe that any other symmetry breaking pattern (such as Z(6) -^ nothing) is 
possible. 

""^^Within a mean field approximation, m^ and m-)^ are the inverse correlation lengths for the real and 
imaginary parts of the trace of the Wilson line, respectively. A ratio of three for m^/m^ at the point 
where the tree-level potential undergoes symmetry restoration is in accord with the mean field analysis 
of Dumitru and Pisarski |2^ for a Z{3>) symmetric model involving only the trace of the Wilson line 
(with up to quartic interactions). 
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on the transition line there will be four co-existing equilibrium states — just as in the 
full SU{3) gauge theory. 

Exploring this two-dimensional parameter space via numerical simulations is surely 
feasible. As a warm-up, one could focus on a single line with a fixed value m'i/v'^, and 
only vary gl/v. It would be very interesting to see if the behavior of physical observables 
(such as the ratios of correlation lengths in different symmetry channels) along such a 
line closely mimics the behavior in hot SU{3) gauge theory as the temperature varies 
from below Tr to far above. 



6. Conclusion 

We have formulated a three-dimensional Z(3)-invariant effective theory which can re- 
produce properties of hot SU{3) Yang-Mills theory both at asymptotically high temper- 
atures, and near the confinement /deconfinement transition. Lowest order perturbative 
matching plus comparison of domain wall properties was used to fix the parameters 
of the effective theory at high temperatures. The form of the effective theory ensures 
that as one increases the dimensionless gauge coupling g'^/v (which increases the size 
of fluctuations), or decreases the heavy mass m^ (which changes the shape of the tree- 
level potential), the theory will undergo a first order phase transition from a deconfined 
phase with spontaneously broken Z{3) symmetry to a confining phase with unbroken 
Z{3) symmetry. This reproduces, by design, the qualitative behavior of hot SU{3) 
Yang- Mills theory near Tc. Numerical simulations of the effective theory will be needed 
to determine the tuning of parameters that makes the effective theory, near its phase 
transition, most closely reproduce quantitative properties of hot Yang-Mills theory 
near the confinement transition. As discussed in the previous section, just varying g'^/v 
while holding the other parameters fixed at the values determined by high temperature 
matching will give a theory that reproduces the correct high temperature behavior at 
small gauge coupling, together with a fluctuation induced phase transition at a critical 
value of gl/v. This, we hope, will be an accurate effective theory for quarkless QCD 
for temperatures ranging from near Tc to asymptotically large. We look forward to 
numerical simulations testing this hypothesis. ^^ 

In the case of a positive result, it will be interesting to investigate generalizations 
of our effective theory to other values of A^c, or to QCD with dynamical quarks. For 
sufficiently heavy quarks, it is known that a weakly first order finite temperature phase 
transition persists, even though fundamental representation quarks break the Z{3) 
center symmetry and no local order parameter for the transition is available. Therefore, 
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Such simulations are currently underway U4 



-24- 



dynamical quarks may be viewed as adding soft Z{3) symmetry breaking terms to the 
effective theory. Such terms will lift the degeneracy of the minima of the potential, 
leaving the minimum with a real expectation value for the trace of the Wilson line as 
the unique global minimum. Matching, at high temperature, to EQCD with a non- 
zero number of quark flavors should be straightforward. The resulting effective theory 
should capture some of the Z[3) physics postulated to be important in the full theory 
even at non-zero N{ [§. Pursuing this generalization will be left for future work. 
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A. Minimizing the potential 

The variation of the hard potential Vo{Z), defined in Eq. ( |3.8| ), is 

6Vo = ti [6Z {ci Z^ + C2 Z~^ det Z + 2c3 Z^ Z Z^}] + h.c. (A.l) 

and this vanishes when Z satisfies 

CiZ^Z + C2detZ + 2c3{Z^Z)'^ = 0. (A.2) 

Subtracting this equation from its Hermitian conjugate shows that det Z must be real 
at extrema. For C2 < 0, which we have assumed, one may see directly from the potential 
(|3.8|) that the minimum of Vo{Z) will occur for matrices with positive determinant. 

Let A^ denote a (positive, real) eigenvalue of the the positive-definite matrix Z'^Z. 
To satisfy relation ( [A .21 ), these eigenvalues must obey 

2c3 A^ + ci A^ + C2 det Z = . (A.3) 

For C2 < < C3, and a positive determinant of Z, this has a unique positive root, 



^2 _ -ci + ^/c^~^^8c^c^detZ 



Acs 

Therefore, all eigenvalues of Z^Z are equal, implying that this matrix is proportional 
to the unit matrix. Consequently, Z must equal a special unitary matrix times A. 
So detZ = A^ and condition ( [OD becomes A^ = [-ci + \/cf^^8c^c^]/{Ac3). The 
solutions of this equation are A = and A = A-t where 



_ -C2 ± y/c^ - 8C1C3 

A± = — . (A.5) 

These non-zero roots are real if c^ > 8C1C3. Evaluating Vq at its three possible extrema, 
one finds that the root A+ corresponds to a global minimum provided C2 < and 
C2 > 9C1C3. Assuming C2 satisfies these conditions, Vo{Z) is minimized when Z = X^Q 
for arbitrary Q G SU{3). When ci > there is a local minimum at Z = and 
the condition cl > 9ciC3 implies that C2 must be less than — Sy^ciCs. But if ci < 
then 2^ = is a local maximum and C2 must merely be negative. If ci > and 
8C1C3 < C2 < 9ciC3, then ^ = is the global minimum of the potential and Z = X^Q 
is only a local minimum. 

To analyze the soft potential Vi{Z), given in Eq. (|3.9| ), it is convenient to insert the 
singular value decomposition M = LmK^ . Here L and R are unitary matrices, while 
m is diagonal, real and positive. Only the cubic terms in Vi{Z) depend on the unitary 
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matrices L and R, and it is easy to see that £2 (tr [M^] + tr [(M"'')^]) is minimized for 
L = R ii C2 < 0, or L = —R if C2 > 0. Therefore, it is sufficient to minimize Vi{Z) for 
M a diagonal real matrix. This gives 

V,{Z) = Y, (gi /x^ + 252 /i? + £3 /i- ) , (A.6) 

i 

where {/ij} are the real eigenvalues of M. Solving for zeros of this expression, one finds 
that Vi{Z) vanishes only at M = 0, and is strictly positive for all non-zero M, provided 
ci and £3 are positive and C1C3 > C2, which we assume. 

Non-trivial minima of the hard potential Vq{Z) [i.e., Z/X^ G SU{3)] coincide with 
minima of the soft potential Vi{Z) [i.e., configurations satisfying 2 oc 1 so that M = 0] 
when Z equals A+ times an element of the center of SU{3). So the global minima of 
the complete potential Vo(Z) + g^ Vi{Z) (when C2 > OciCs) lie at 

Z = e2™/3A+1, ne {0,1,2}. (A.7) 



Using the decomposition ( p.l2| ) of Z and the resulting relations (|3.15| ) between 
the coefficients in Vo{Z) and the tree level heavy masses, one may translate the above 
conditions on the coefficients into equivalent conditions on m^ and m^ (which are the 
tree level masses for fluctuations in the magnitude and phase of tiZ, respectively). 
The results may be summarized as follows. 

1. If 171"^ > and m^ + 3r7i| > then C2 < < C3. These conditions are required so 
that Vq{Z) is bounded below, and favors matrices with positive determinant. 

2. If m^ > wi?,/3 then ci is negative and there is a local maximum ai Z = and 
three degenerate global minima at Z = \^ e^'^*"'^l. 

3. If m^/9 < m| < rn'^/3 then ci > and C2 < —Sy/cic^, implying that there is a 
local minimum at Z = and three degenerate global minima at Z = \^ ^2mn/?,^_ 

4. If < m^ < m^/9 then ci > and — 3^CiC3 < C2 < — v^SciC^, in which case 
there is a global minimum at ^ = and three degenerate local minima at ^ = 

B. Expansion in shifted fields 



Expressing the potential (|3.71 ) in terms of the shifted flelds ( |3.12|) which describe fluc- 
tuations away from the minimum gives 

gf V{Z) = Kni, + \ml <\? + \ml x' + ml tx{h^) + V,^,{<\>, x, h, a) , (B.l) 

with 
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K.in = ^ (9ci + VC2) v^/gl , (B.2) 

and 
VinM, X, h, a) = -^ I (3ci + hc2) (p^ + (3ci + 2vc2) 4>x^ + (6ci + vc2)4>tT (a^) 

+ (18ci + 7vc2) (/)tr (/i^) + (12ci + 2i;c2) xtr (a/i) 
+ Vg (6ci + 4f C2) tr (a^h) - Vg (6ci + ^VC2) tr (/i^)! 

- % (3ci + t;c2) { I (0' + x') ' + |0' tr (a^ + 3/1^) + 20x tr (ah) 

+ \x^ tr (3a^ + /i^) + 20 tr [c?h + /i^) + 2x tr (a^ + a/i^) 
+ I tr (a^ + 4a2/i2 _ ^{aKf + /i^) } 

+ ^gjci tr (a^ + /i^)} - ^|{2c2tr (30^/1 - h^)} 

+ (7^{c3tr(a^ + 4a2/,2_2(a/i)2 + /i4)}, (B.3) 



Here C3 has been expressed in terms of v using Eq. (|3.11|) . The tree level masses are 



m| = -2ci - |f C2 , rr?^ = -vc2, m^ = -2ci - |t; C2 . (B.4) 

Expressed in terms of these masses, \4iin = ~ los (^"^0 ~ "^-x)! 9\- This shows directly 
that rr?, > ^rn^ is necessary for the global minima to lie at the non-trivial Z{3) extrema. 
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